Subcritical temperature in Bose - Einstein condensates of chiral molecules 
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Bose-Einstein condensation of a sample of non-interacting chiral molecules leads to a non zero 
optical activity of the condensate and also to a subcritical temperature in the heat capacity. This is 
due to the internal structure of the molecule which, in our model, is considered as a simple two-state 
system, characterized by tunneling and parity violation. The predicted singular behavior found for 
the specific heat, below the condensation temperature, sheds some light on the existence of the so 
far elusive parity violation energy difference between enantiomers. 
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Since the prediction and subsequent discovery of parity 
violation [H in weak interactions, the role of discrete 
symmetry breaking in fundamental physics is an intrigu- 
ing field of research. Although the weak interaction (be- 
tween electrons and nuclei) mediated by the gauge Z° bo- 
son has been extensively studied and observed in atoms 
Q, it has only been predicted in molecules. The impor- 
tance of this symmetry-breaking in molecules is twofold: 
(i) it could be intimately related to the origin of ho- 
mochirality, that is, the almost exclusive one-handedness 
of the chiral molecules found in living systems, this be- 
ing one of the most fascinating open problems which 
links fundamental physics with the biochemistry of life 
[3] and (ii) at fundamental level, the intrinsic chiral na- 
ture which is present in some molecules should reflect 
the underlying interaction containing pseudoscalar mag- 
nitudes as those appearing in the weak one. In molecular 
systems, the theory of electroweak interactions predicts 
a parity violating energy difference (PVED) between the 
two enantiomers of chiral molecules to be between 10"^'^ 
and 10^^^ eV (5l-[TT|. However, no conclusive energy dif- 
ference has been reported, for example, in experimental 
spectroscopic studies of the CHBrClF molecule reaching 
an energy resolution of about 10~^^ eV [l^ [isj. 

Due to the fact that the parity violating signals are 
easily masked by thermal effects, it is highly desirable 
to reach cold or ultracold regimes in the laboratory [l^ - 
[l7j . The measurement of fundamental physical proper- 
ties such as the PVED between enantiomers could be 
achieved by trapping molecules at low temperatures (in 
the milikelvin range or below) and, subsequently, per- 
forming ultrahigh-resolution spectroscopic measurements 
of vibrational or electronic transitions. We would like to 
point out that these temperatures, and even colder, have 
been reached within the actual experimental capabilities 
only for atoms and diatomic molecules. The rapid ex- 
pansion of the research field of ultracold chemistry is al- 
ready opening new and exciting possibilities concerning 
more complex systems [l8| . Thus, the study of quantum 
thermal effects in chiral molecules through pseudoscalar 
operators as the basic object describing their thermody- 
namics is of fundamental importance. Very recently, we 



have studied the corresponding classical thermodynam- 
ics of non-interacting chiral molecules [l^, HSl- In this 
sense, as pointed out by Flambaum and coworkers con- 
cerning some capabilities of Bose-Einstein condensates 
to amplify weak interactions in chiral molecules [l^ , the 
study of Bose-Einstein condensation (BEG) for a gas of 
chiral molecules could shed some light on some properties 
relevant to determine the so far elusive PVED. 

On the other hand, the study of cold atoms with in- 
trinsic spin structure is being a very hot topic, includ- 
ing two-component condensates, spinor condensates, etc 
|2l| . In this sense, for a Bose gas of chiral molecules, the 
inclusion of extra information about the internal struc- 
ture of such molecules could be determinant to obtain 
some macroscopic properties relevant to determine the 
PVED such as, for example, by measuring the optical 
activity and heat capacity. For this end, we will employ 
a two-state model which includes tunneling and parity 
violating effects in chiral molecules (this approach was 
first considered in and is widely used to describe, for 
instance, the internal dynamics of chiral molecules). Al- 
though the two-state (or spin 1/2) model is well known 
in many fields, including spin tunneling, magnetization, 
quantum dissipative systems, to name only a few, this 
is the first time, to the best of our knowledge, that this 
model is employed to study BEG of chiral molecules. Fol- 
lowing the seminal works on BEG, predicted for a non- 
interacting gas and carried out in numerous dilute atomic 
23| and in some diatomic molecules as in 
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gases 

we will assume that the interaction among molecules is 
negligible at first order and, therefore, a non-interacting 
quantum gas of chiral molecules will be considered. 

In this Letter we show that a condensed gas of non- 
interacting chiral molecules displays optical activity. In 
addition, we show a dramatic change of the heat capac- 
ity with temperature, leading to the appearance of a sec- 
ondary maximum (developing a shoulder-type-structure) 
associated with the Schottky anomaly below the conden- 
sation temperature. This subcritical temperature is ex- 
pected to lie in the cold or ultracold regime. 

In an ideal Bose-Einstein gas, only the translational 
part of the molecular motion is involved and the criti- 
cal temperature for BEG is reached when the fugacity of 
the gas z = exp(/3/^), with fi being the chemical poten- 
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tial and /3 = {kBT)~^ , is approaching unity (assuming; 
that the energy of the ground state is zero) [2^- This 
Tc depends on the molecule mass and the density of the 
system. At T < there must be a non- negligible frac- 
tion of molecules in the ground state, the internal energy 
goes with T^/^ and the heat capacity with T^/^. If each 
chiral molecule in the condensed phase has its own inter- 
nal structure, we could ask ourselves if the correspond- 
ing quantum statistics describing the thermodynamical 
behavior, intimately related to the PVED, can be mani- 
fested by a macroscopic effect. For this goal, a two-state 
model is assumed for each chiral molecule. 

In the absence of parity violation, the dynamics of a 
chiral molecule is well described by a symmetric double 
well potential. In this case, the true stationary states 
of a chiral molecule arc the (achiral) eigenstates of par- 
ity. The effect of including in the Hamiltonian the in- 
ternal P-odd term given by the parity violating electron- 
nucleon electroweak interaction, H^^ , introduces an en- 
ergy difference between the two minima of the double 
well potential, leading to a new set of energy eigen- 
states. Let us consider the total Hamiltonian of the sys- 
tem H = + H^^ , with including only parity con- 
serving terms. Using a two-state model, wc choose the 
left and right chiral states basis, \L) and \R) (localized 
respectively in the left and right minimum of the dou- 
ble well potential), in order to show clearly the parity 
properties of the Hamiltonian (2^ : 

epvCTz (1) 



where ax,z are the Pauli matrices. The eigenstates of 
H, 1 1,2), can be expressed as linear combinations of 
the chiral states by means of a rotation where the cor- 
responding mixing angle, 6, obtained from the knowl- 
edge of the eigenvectors oi H , is given by tan 26* = 
{2Hlr/Hrr - Hll) = {5/epv). The energy splitting 
between the two eigenstates of is 2i5 > 0, where 
i5 = {L\H^\R). This magnitude is related to the height 
of the barrier of the double well potential and is in- 
versely proportional to the tunneling time. The PVED 
is then given by \Hrji — Hll\ ~ |2epy| with tpy = 
{R\HP^\R) = ~{L\HP^\L) (we remark that epy is the 
eigenvalue of a pscudoscalar operator in the chiral basis). 
The eigenvalues of the system are given by Ei^2 — -E^oT A, 
with Eo = {Hrr + Hll)/"^ and A = ^Cpy + ^"^ (here- 
after we will take £^o = for the sake of simplicity). 

Bose-Einstein condensation of a non-kinetic-two-state 
gas. In a gas of non interacting chiral molecules at tem- 
perature T, the inclusion of P-odd effects leads to deal 
with a pscudoscalar operator for obtaining any thermo- 
dynamic variable [lol [20l|. Thus, for a system described 
by a biased double well potential, the only magnitude 
which, roughly speaking, distinguishes between left and 
right conformations is the population difference between 
both wells, N]^ — Np, which is directly related to the 
optical activity of the system. If X is a pscudoscalar 
operator, it can be shown that 



where ±x are the eigenvalues of X and the subscript 
stat stands for the different statistics considered. Eq. 
([2]) is a generalization to any statistics of the Maxwell- 
Bolztmann (MB) thermal average [l^. A straightfor- 
ward calculation gives the expression for N2 — Ni to be 
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where p = applies for MB, p = —1 for Bose-Einstein 
(BE) and p = +1 for Fcrmi-Dirac statistics. We remark 
that, when epy dominates over 6, the true eigenstates of 
the system are the chiral states, \L) and \R), and chang- 
ing the basis from |1), |2) to \L), \R) introduces the factor 
cos 261 = epv/A leading to Nl-Nr = epv(7V2 - 7Vi)/A. 
Moreover, — Nr changes its sign when considering 
a parity-transformed double well since epy changes to 
—epy. From the knowledge of N2 — Ni, the internal en- 
ergy of the gas is 



U = -A{N, - N2)s 



(4) 



From Eq. ^ one can inmediately obtain the heat ca- 
pacity at constant volume. This Gu can be interpreted 
as a measure of the fluctuations of the optical activity of 
the system, that is, a measure of the fluctuations of the 
pscudoscalar character of the system. 

Let us now consider BE statistics (p = — 1) for a non- 
kinetic-two-state gas, that is, the Hamiltonian H includes 
tunneling and parity violation but not kinetic terms. The 
total number of particles is N = Ni + N2, where Ni^2 is 
the number of particles in the ground, |1), and the excited 
state, 1 2), respectively. Using the standard expression 
for the occupation numbers for Bose statistics, we get 
N = Ni-\- (z-iexp(/3A) - 1)-^. At the condensation 
temperature, T*, all the particles are in the excited state, 
N2 ■ (this definition of the condensation temperature was 
also pointed out in [1^). This leads to the condensation 
temperature to be determined from 
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The condensate displays a non-zero optical activity 
which is proportional to 



N2- Ni = l 



exp(2;3*A) - 1 



exp(2/3A) - 1 ' 

and this enables us to write, for T < T*: 

No ^ exp(2/3*A) - 1 
N exp(2^A) - 1 ' 

the average internal energy as 

U cxp(2^*A^-2 
" ^ ]V ^ exp(2/3A) - 1 - ^) 

and the heat capacity as 
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It is worth noting that this heat capacity reaches its max- 
imum value at the subcritical temperature, Tg^, 



l3scA ~ 0.797 



(10) 



which reflects nothing but the Schottky anomaly due to 
the saturation of the energy levels of the system. 

Bose-Einstein condensation of a kinetic-two- state gas. 
If we include a kinetic term in the Hamiltonian, the total 
partition function for an ideal Bose gas of chiral molecules 
in the two-state model can be written as 



(11) 



where the subscripts denote the kinetic and internal con- 
tributions, respectively. The first factor can be found in 
any standard textbools of Statistical Mechanics [1^ and 
the second factor is given, for our case, by 
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The factorization of the partition function leads to a sum 
of the corresponding heat capacities, C*°* = C^*" -|-C*"*. 
We note that there are now two temperature regimes in- 
volved in the chiral system: (i) the temperature at which 
the gas undergoes BEC (Tc) and (ii) the subcritical tem- 
perature (Tsc, below Tc) where the heat capacity displays 
a maximum. In the case of a free gas, Tc is related to the 
mass and number density of the bosons but T^c depends 
only (under the two-state model) on the energy splitting 
between |1) and |2) states. Thus, for T < T^, C*"* is 
given by Eq. ^ and, for T > Tc, the corresponding 
expression can be obtained from Eq. (g]). C^™ is the 
usual heat capacity for an ideal Bose gas. It is worth 
pointing out that the very interesting case where Tc and 
Tsc are of the same order gives place to an appreciable 
change in the heat capacity and to a shoulder-type struc- 
ture of it below Tc- This behavior is displayed in Fig. [1] 
where the total heat capacity is plotted in terms of the 
reduced temperature ksT/A for the cases Tsc = O.lTc 
(top panel) and Tsc = 0.5Tc (bottom panel) (we have 
assumed T* = Tc). 

Strictly speaking, the effect of A on the condensation 
temperature, T^ results into a very small shift. Thus 
for a kinetic-two-state-gas we get that the corrected or 
shifted condensation temperature is given by 
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This difference is so tiny that for most practical purposes, 
we keep Tc instead of T. 

When Tc = Tsc, the subcritical temperature corre- 
sponds to a very precise density of bosons of mass m. 
Use, given by 
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Figure 1. Total heat capacities for a gas of chiral molecules 
as a function of the reduced temperature for two cases: Tac = 

0. 1Tc and Tsc = O.STc. The gas undergoes BEC at fcsT/A = 

1. We have taken T* — Tc. See text for symbols. 



The bottom panel of Fig. ^ should be observed, for 
example, for the T2Se2 system. This system has epy ~ 
2.5T0-^* eV, i5 « STO"^'' eV and the tunneling time is 
~ 40 s (see Table 2 of [IH). Its subcritical density, n^c, 
is about 10^" cm-^^ which is of the order of the critical 
densities. Thus, in this case, the corresponding disconti- 
nuity in the heat capacity could provide us a direct and 
clear signal of molecular parity violation. 

In Fig. (12); a three dimensional plot of the heat ca- 
pacity in the range of variation [0,1.2] of the reduced 
temperature and [10^^^, 10^^^] of the splitting A (eV) is 
showed. In this plot it is clearly seen the overall shoulder- 
type structure as well as the subcritical maximum for dif- 
ferent chiral molecules, characterized by the splitting A. 
This splitting is representative of molecules such as, for 
example, H2S2, CHBrClF and H2Se2 [ll|. 

The subcritical maximum (or anomalous shoulder-type 
structure) found here could be considered as a signal of 
the Schottky anomaly which comes from lifting the de- 
generacy of the internal degrees of freedom. Although 
this was previously considered in the framework of a 
path integral approach to study BEC in spinor gases 
[27|, our simple model allow us to identify more clearly 
the anomaly. Wc would like to stress that the subcritical 
temperature is determined in essence by 6 and epv, these 
temperatures lying in the cold or ultracold regime. The 
existence of that temperature could be possible for any 
system whose constituents have internal ground states 
slightly splitted, such as for example: the inversion dou- 
bling of non-rigid pyramidal molecules, the splitting due 
to torsional or internal rotations through potential bar- 
riers in non-rigid molecules, or eventually the hyperfine 
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Figure 2. Three dimensional plot of the heat capacity versus 
A and the reduced temperature. The gas undergoes BEC at 
ksT/A = 1. We have taken T* — Tc- See text for symbols. 



structure of the ground state of atoms and molecules (in 
this last case, it could be more easy to verify the type 
of anomaly here studied). However, chiral molecules in 
this two-state model play a special role in BEC when the 
splitting is mainly due to parity violation (as in T2SC2). 
In this case, an observation of the anomaly predicted in 
the heat capacity provides a direct information about 
the PVED. In addition, when epv determines the dy- 
namics, then the eigenstates of the system tend to be the 
chiral states \L,R). We also point out that a measure- 
ment of the optical activity of the condensate would be a 



confirmation of the existence of the PVED, noting that 
the factor epy/A, which is the optical activity for zero 
temperature, should be measurable by state-of-the-art of 
actual polarimeters, as we concluded in [20j . 

Finally, we would like to point out that the energy scale 
of molecular parity violation is associated with the natu- 
ral scale of densities which corresponds to those achieved 
in Bose-Einstein condensates. Thus, BEC seems to be 
essential to detect the PVED. Furthermore, as is it well 
known, interactions could modify the physics of the prob- 
lem so it is fundamental to ask ourselves to what extent 
the anomaly in the heat capacity persists when we con- 
sider an interacting system of chiral molecules. Although 
one could introduce the interaction between molecules 
under the Gross-Pitaevskii or any other more sophisti- 
cated treatment, an alternative way is to study the molec- 
ular sample as an open quantum system under the influ- 
ence of dissipation. In this sense, we note that the heat 
capacity anomalies of open quantum systems have been 
recently studied when coupled to a thermal bath j28j . 
showing the robustness of these anomalies when the in- 
teraction is taken into account. Thus, we expect that 
a similar stability for the anomaly here considered will 
persist even in presence of an environment. The study 
of heat capacity anomalies in BEC chiral gases embed- 
ded in appropiately chosen environments is currently in 
progress. 
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